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Abstract 

We present an estimate of the cross-section for the exclusive production of a 
p^-meson pair in e + e~ scattering, which will be studied in the future high en- 
ergy International Linear Collider. For this aim, we complete calculations of the 
Born order approximation of the amplitudes Yl t(QT)1l riQl) ~~ * P°lP°li f° r arbi- 
trary polarization of virtual photons and longitudinally polarized mesons, in the 
kinematical region s ^> —t, Q\ , Q\. These processes are completely calculable in 
the hard region Q\ , Q\ ^> Aq CD and we perform most of the calculations in an 
analytical way. The resulting cross-section turns out to be large enough for this 
process to be measurable with foreseen luminosity and energy, for Q\ and Q\ in 
the range of a few GeV 2 . 
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1 Introduction 



The next generation of e + e~— colliders will offer a possibility of clean testing of QCD 
dynamics. By selecting events in which two vector mesons are produced with large ra- 
pidity gap, through scattering of two highly virtual photons, one is getting access to the 
kinematical regime in which the perturbative approach is justified. If additionally one 
selects the events with comparable photon virtualities, the perturbative Regge dynamics 
of QCD of the BFKL [I] type should dominate with respect to the conventional partonic 
evolution of DGLAP [2] type. Several studies of BFKL dynamics have been performed 
at the level of the total cross-section [31 H] . 

In the paper [5] the diffract ive production of two J/^> mesons was studied as a promis- 
ing probe of the BFKL effects. Recently, we have advocated P, [3, [8] that the electro- 
production of two p— mesons in the 7*7* offers the same advantages. In this case the 
virtualities of the scattered photons play the role of the hard scales. A first step in this 
direction was made by considering this process with longitudinally polarized photons and 
p— mesons, 

il(qi) i*M - pl(ki) P l(k2) , (1.1) 

for arbitrary values of t = (q\ — ki) 2 , with s ^> —t. The choice of longitudinal polarizations 
of both the scattered photons and produced vector mesons was dictated by the fact that 
this configuration of the lowest twist-2 gives the dominant contribution in the powers of 
the hard scale Q 2 , when Q\ ~ Q\ ~ Q 2 . 

The aim of this study is to complete the Born order evaluation of the cross-section of 
the process 

e + e~ ^e + e-p L p° L , (1.2) 
illustrated in Fig{TJ We calculated contributions with remaining combinations of polar- 
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Figure 1: Amplitude for the process e + e — > e + e p° L p° L . 
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izations of virtual photons necessary to obtain all helicity amplitudes of the processes 

t£,t(9i) - fith) P ° L (k 2 ) . (i.3) 

Let us note, that the double tagging of final leptons gives in particular the possibility to 
separate the contributions of various photon polarizations, entering in ( 11. 2\\ and thus to 
study the corresponding parts of the cross-sections which are computed in this paper. We 
are focusing here on the high-energy limit in which t— channel gluonic exchanges dominate. 
On the other hand, in the description of the process (ll.2p . there is a potential possibility 
that also contributions with rather small s 7 * 7 * have to be taken into account. In this 
case one should include in principle both quark and gluon exchanges. The contribution 
of quark exchange was analyzed in [9]. This quark-box contribution is investigated in 
subsection 14.31 

We will also not consider here the case of transversally polarized p— mesons. It would 
require to deal with possible breaking of QCD factorization [10J, although a method to 
overcome this problem has been proposed [TT] . 

The BFKL enhancement was studied for t — in [5] and [12]. In this latter case, 
the peculiar value t — automatically selects the longitudinally polarized photon. A 
dedicate study for arbitrary value of t should thus be performed to get an evaluation of 
BFKL enhancement effects of the Born order evaluation performed in the present paper 
for transversally polarized photon. This problem will not be addressed here. 



2 Kinematics 

The measurable cross section for the process (11.21) of Fig{T] is related to the amplitude 
of the process (jl.3p . illustrated in FigfJl through the usual flux factors for respectively 
transversally and longitudinally polarized photons 

Kvi)= l + i(w) = i-%, (2.1) 

where yi (i = 1, 2) are the longitudinal momentum fractions of the bremsstrahlung photons 
with respect to the incoming leptons. This relation reads [T3] 

da (e+e~ -> e + e-p° L p° L ) 
dy x dy 2 dQ\dQ\ 

= ^2 (-)' [Kvi) Ifa) - P°lP°l) + t(yi) l(y 2 ) a( 7 T7l - P°lP°l) 

+ l( Vl ) t(y 2 ) - pIpD + t( yi ) t(y 2 ) <x( 7 ^7t - P°lP°l)} ■ 

The presence of hard scales Q\ permits us to apply the collinear approximation at each 
qqp— meson vertex, and the use of distribution amplitude (DA) for describing the qq 
content of the p mesons, as illustrated in FigJ21 In this paper, except for section 14.31 the 
amplitude Mh will be described using the impact representation, valid at high energy, as 
illustrated in Figj3l 
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Figure 2: The amplitude of the process 7*(Qi)7*(Q2) — ► Pl(^)p°l(.^) with the collinear fac- 
torization in the qqp vertex. 



Let us introduce two light-like Sudakov vectors q[ and q 2 which form a natural basis 
for two scattered virtual photons, which satisfy 2q[ ■ q' 2 = s ~ 2qi • q 2 - The usual s 7 * 7 * is 
related to the auxiliary useful variable s by s 7 * 7 * = s — Q\ — Q\. The momentum transfer 
in the t— channel is r = k± — q±. In this basis, the incoming photon momenta read 

/ Ql 1 1 I Q'2 1 /r, o\ 

?1 = ?1 ?2 alld ?2 = ?a ?1 • ( 2 - 3 ) 

s s 
The polarization vectors of longitudinally polarized photons are 

^ (1) = ^ + ^4 and e f) = ^ + ^, (2.4) 

with e|™ = 1 and q^ ■ = 0, whereas the polarization vectors of transversally polarized 
photons are two dimensional transverse vectors satisfying — —1 (i — 1, 2) and qi ■ 
£r(i) = 0. 

We label the momentum of the quarks and antiquarks entering the meson wave func- 
tions as l\ and l[ for the upper part of the diagram and I2 and l 2 for the lower part (see 
FigEj). 

In the basis ( 12. 3p . the vector meson momenta can be expanded in the form 

rr .2 



kx = a(ki) q[ + q 2 + r± 

ac(ki) s 



k 2 = f3{k 2 )q' 2 + —f—q' 1 -T ± . (2.5) 
p{k 2 ) s 

Note that our convention is such that for any tranverse vector v± in Minkowski space, v 
denotes its euclidean form. In the following, we will treat the p meson as being massless. 
a and (3 are very close to unity (explicit expressions can be found in [7j), and reads 

s \s z 

(3(k 2 )^l-9l±^ + (\) ^ (2.6) 
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Figure 3: The amplitude Mh in the impact representation. The vertical blob symbolizes the 
interaction of two qq dipoles through gluon exchanges at high s. 

where r 2 = — r\. They will be replaced by 1 in the phenomeno logical applications of 
sections H] and [5j In this decomposition, it is straightforward to relate t = r 2 to r 2 
through the approximate relation 



(see (7j for an exact relation). From Eq. (12.71) the threshold for \t\ is given by \t\ m m — 
QIQ2/S, corresponding to r_|_ = 0. In the kinematical range we are interested in, the 
relation (12. 7p can be approximated as r 2 = —t, as usually in the Regge limit. 

The links with the e + e~ process can be made by using the same Sudakov basis for the 
two incoming leptons: 

pi = —q'i +yi—Q2 +p±i and P2 = —q' 2 +V2—q'i +P±2 , with p 2 = — ^Q 2 . (2.8) 

yi s y 2 s yf 

Thus, one gets 



In the rest of the paper, since we keep only the dominant s contribution, we use the 
approximate relation s e + e - ~ s/iy^)- 

3 Impact representation 

The impact factor representation of the scattering amplitude for the process (11.31) has the 
form (see FigH]) 




(2.7) 
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Figure 4: Amplitude Mjj at Born order. The t-channel gluons are attached to the quark lines 
in all possible ways. 

/72 7 
(27r) 4 fc (r — A;) 2 

where J^ 9l) ^ kl \k,r- k) (jT^ (92Hp ° (fca) (fc, r - /c)) are the impact factors corre- 
sponding to the transition of 7£t(<?i) — * Pl(^i) (7lt(<?2) — ► PzX^)) via the t— channel 
exchange of two gluons. The amplitude (13.11) calculated in Born order depends linearly on 
s (or s 7 * 7 * when neglecting terms of order Qj/s) as the impact factors are s- independent. 

Calculations of the impact factors in the Born approximatiorQare standard [H]. They 
are obtained by assuming the collinear approximation at each qqp— meson vertex. Pro- 
jecting the (anti)quark momenta on the Sudakov basis q[, q' 2 , 

h = Zxq[ + Ixi + zxrx. 
l[ = ziq[ - Ixi + zxrx 
h = z 2 q' 2 + l±2 ~ z 2 r±_ 
1' 2 = z 2 q' 2 - I ±2 - z 2 r±_ 

we put the relative momentum Zj_|_ to zero. For longitudinally polarized photons the 

1 Now, the forward impact factor of 7£(Q 2 ) — > p\ transition has been calculated at the next-to-leading 
order accuracy in [To] . 



(J±l + Zir±) 2 , 
Q 2 > 

Z\S 

(-Zj.i + z x r±) 2 , 
Q 2 ) 

Z\S 

(I ±2 ~ z 2 r ± ) 2 , 

Qi , 

z 2 s 

t±^Z™±£ q > , (3.2) 

Z 2 S 
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impact factor reads 



JTlM^Mfar-k) = S^a^^Qifpaiki) J dziZi z { (f>(zi)P P (z h k, r, m) , (3.3) 

o 

where the expression 

P P (zi,k,r, Mi ) = z 2 l2 + /i 2 + z2 £2 + // 2 ~ ( Zi r - k) 2 + /i 2 ~ (zir-k) 2 +//? (3 ' 4) 

originates from the impact factor of quark pair production from a longitudinally polarized 
photon. 

For transversally polarized photons, one obtains 

i 

e 5 ab f 

= Aii 2 a s -^= — f p a(ki) J dzi (zi - z { ) cp(zi) e ■ Q{zi, k, r, fa) , (3.5) 



where 

— 3 fr + nf zfr 2 + fif [ZiZ - k) 2 + \i\ [Zi r-k) 2 + fi- 

is proportional to the impact factor of quark pair production from a transversally polarized 
photon. 

In the formulae (13. 4p and (I3.6P and for the rest of the paper, we denote pf = Q 2 Zi z~i + 
m 2 , where m is the quark mass. The limit m —>■ is regular and we will restrict ourselves 
to the light quark case, taking thus m — 0. Both impact factor (13 .4p and (13.61) vanish 
when fc-»0orr-fc^0 due to QCD gauge invariance. 

In the formulae (13.31 13.5j) . 4> is the distribution amplitude of the produced longitudi- 
nally polarized p°— mesons. For the case with quark q of one flavour it is defined (see, e.g. 
|16j ) by the matrix element of the non-local, gauge invariant correlator of quark fields on 
the light-cone 



(0\q(x) Y l q{-x)\ PL {p) = qq) = f P j dz e ^-D0") , (3.7) 

o 

where the coupling constant is f p — 216 MeV and where the gauge links are omitted to 
simplify the notation. <fr is normalized to unity. The amplitudes for production of p°'s are 
obtained by noting that \p°) = l/\^2(\uu) — \dd)). 

Note that Eq. ( 13.71) corresponds to the leading twist collinear distribution amplitude. 
Such an object can be used strictly speaking for asymptotically large Q 2 . In the phenome- 
logical application of sections 14.21 and 15.31 in order to get measurable cross-sections, the 
dramatic decrease of the amplitudes with increase of Q 2 , combined with the experimental 
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condition of ILC project, requires rather low values of Q\ (of the order of 1 GeV 2 ) for 
which subleading twist contributions could be significant. This can be taken into account 
within a more phenomenological approach which incorporates intrisic kx quark distribu- 
tion and which goes beyond standard QCD collinear factorization [17J. In the present 
paper we do not consider these effects and adhere to the collinear QCD factorization. 

Let us label the amplitudes for the scattering process (11.31) through the polarization 
of the incoming virtual photons as -MaiA 2 - They can be calculated using Eqs.1 13.ip and 
Eqs. (I3.3tl3.6j) supplemented by the choice of the transverse polarization vectors of the 
photons 

e± = -^=( T 1, -i) (3.8) 

and the longitudinal polarization vectors (12 Ah . For the case Ai = A 2 = : 

i 

M o = isCQiQ 2 dz 1 dz 2 z 1 z 1 (j)(z 1 ) z 2 z 2 (f)(z 2 )M 00 (z ll z 2 ) , (3.9) 



with 

/d 2 k 
k 2 ( r Ik) 2 Pp ( Zl ' -' -' ^ Pp ( Z2 ' ~-' ~-' ; ( 3 ' 10 ) 



for the case A2 = +, — : 

1 

C f 

M \ 2 =is — Q l dz 1 dz 2 z 1 z 1 <f){z 1 )(z 2 -z 2 )<f)(z 2 )M QX2 (z 1 ,z 2 ), (3.11) 


with 

/d 2 k 
k 2, ~ k y p p( z i> r, fix) Q(z 2 , -k, -r, /i 2 ) ■ e x ' 2 ; (3.12) 

for the case Ai = +, — : 

1 

C f 

M Xl o = is — Q 2 ldz 1 dz 2 (zi-zi)<p(zi)z 2 z 2 (f)(z 2 )Mx 1 o(z l ,z 2 ), (3.13) 



with 

/d 2 k 
k 2^ r Z k y Q( Zl ' k ^ ^1) • i Al Pp(z 2 , -k, -r, /i 2 ) . (3.14) 

and for the case Ai = +, — , A2 = +, — : 

1 

C f 

M Xl x 2 =is— ldz 1 dz 2 (z 1 -z 1 )<p(z 1 )(z 2 -z 2 )<p(z 2 )Mx 1 \ 2 (z 1 ,z 2 ), (3.15) 


with 

M AlAa (*i, z 2 ) = J k2( f^ k)2 Q(zi, k, r, n x ) ■ e Al « Q(z 2 , -k, -r, /i 2 ) ■ e x ^ . (3.16) 
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Here and in the rest of this paper, we denote C = 2 n a 2 s a em / 2 . In terms of the 
above amplitudes, the corresponding differential cross- sect ions can be expressed in the 
large s limit (neglecting terms of order Qf/s) as 



tfc/r^-pM \M Xl x 2 \ 



dt 16 7T s 2 

and it does not depend on s. 



(3.17) 



4 Non-forward Born order differential cross-section 

for 7i,r 7i, T Pl Pl 

4.1 Analytical results for /c^-integrated amplitude M\ ± \ 2 

In this section we summarize the results for the amplitudes M\ 1 \ 2 obtained after perform- 
ing analytically the k± integrals. Such analytic expressions give us the effective possibility 
of studying various kinematical limits in the variables Q\ , Q 2 ., t. The k± integrations were 
done using the method of Ref.fTj which exploits in an efficient way the scaling proper- 
ties of integrals appearing in conformal field theories. The generic k± integral involves 
an integrand which corresponds to a box diagram with two distinct massive propagators 
and two massless propagators. Because of that, the k± integrations result in long and 
complicate expressions. Thus, we discuss below only the general structure of the results 
and we relegate all technical details of k± integrations to the Appendix. 

In the transverse-transverse (TT) case, the amplitude can be expressed in term of two 
projection operators in the transverse plane as follows: 



a(r, Qi, Q2; zi, z 2 ) 6 lJ — + b(r; Q 1? Q 2 ; z 1: z 2 



e- 1 ^, (4.1) 



where we denote r 2 = r 2 . 

Combining (13.151) and fl4.ll) . and using |M ++ | 2 = |M | 2 , one gets in the case of two 

photons with the same polarization : 

\M ++ f = \M^\ 2 = s^ (4.2) 



x 



1 

dz x dz 2 (21 - z x ) <t>{zx) {z 2 - z 2 ) (p(z 2 ) (6(r; Qx, Q 2 ] Z\, z 2 ) - a(r; Q u Q 2 ; zx, z 2 )) 







and analogously for different polarizations : 

l-M + -| 2 =* 2 ^ (4.3) 



x 



1 

/ dz x dz 2 (zx - Zx) (j){zx) (z 2 - z 2 ) <p{z 2 ) (6(r; Qx, Q2; z x , z 2 ) + a(r; Qx, Q2; Zx, z 2 )) 
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For the longitudinal-transverse (LT) case, restoring the dependency over all variables, one 
defines from (13.121) and (13.141) the scalar function / 

M A(r; Qi, Q2; z u z 2 ) = /(r; Qx, Q 2 ; z 1} z 2 ) r ■ e A , (4.4) 

or equivalent ly 

M A ofc Qi, Q2; Zi, z 2 ) = /(r; Q 2 , Qx; z 2 , z x ) r • e A , (4.5) 

which leads to 



\M 0+ \ 2 =\M ^ = s 2 ^QlrJ 



dzi dz 2 z\ zi (f)(zx) (z 2 - z 2 ) (f>(z 2 )f(r; Qx, Q 2 ; Zx, z 2 ) 

(4.6) 



and analogously for the transverse-longitudinal (TL) case 



\M +0 \ 2 = \M- \ 2 = s 2 ^QjrJ 



dzi dz 2 z 2 z 2 4>(z 2 ) (zx - zx) (f>(zx)f(r, Q2, Qi] z 2 , z x ) 



(4.7) 

The expressions of a(r; Qx, Q 2 ; z 1: z 2 ), b(r;Qx,Q 2 ; z 1: z 2 ) and /(r; Qi, Q 2 ; £i, z 2 ) pre- 
sented as combinations of finite standard integrals are given in the Appendix. 

For the longitudinal- longitudinal (LL) case, it turned out [7] that ( 13.91) can be effec- 
tively replaced by M{z x , z 2 ) whose integral over zx, 2 with symmetrical DA gives the same 
result. M(zx, z 2 ) reads 

Mqq{Zx, Z 2 ) = - ( 2 2 2 + _ 2 2 i 2 J ^2(^2) ~ ( 2 2 1 2 + -2 2 1 2 J ^1(^1) 

\ z{r_ l + p\ zfrf + n\ J \ z^r 2 + p 2 z 2 rf + \x 2 J 

) • (4-8) 

J 3m and J4 MlM2 are two dimensional integrals with respectively 3 propagators (1 massive) 
and 4 propagators (2 massive, with different masses), they are both IR and UV finite. 
Their expressions are given in the Appendix. 

Due to the collinear conformal subgroup SL(2,R) invariance [TS], the pi distribution 
amplitude has an expansion in terms of Gegenbauer polynomials of even order which reads 

oo 

4>{z) = 6z(l - z) (1 + J2 ^nC 3 2 i 2 (2z - 1)) . (4.9) 

n=l 



Except for a short discussion in section I5.3[ we restrict ourselves to the asymptotical 
distribution amplitude corresponding to a 2n = 0. 

To complete the evaluation of the amplitude A4, one needs to integrate over the quark 
momentum fractions Z\ and z 2 in the p mesons. For arbitrary values of t, it seems not 
possible to perform the Z\ and z 2 integrations analytically. We thus do them numerically. 
We observe the absence of end-point singularity when zx( 2 ) — > or Zxm — > 1. Indeed, 
for the longitudinal polarizations involving Pp as defined in Eq. (l3.4p . the z divergency of 
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type 1/z, 1/(1 — z) is compensated by the zz factor when z — > 0, 1, while for transverse 
polarizations, involving Q as defined in Eq. (13.61) . there is no singularity since Q is itself 
regular. 

For the special case t = t m i n (where only the LL amplitude is no n- vanishing) , which 
will be useful in the discussion of sections 14.21 and 15. 3[ the integration over Zi can be 
performed analytically, with the result [7] 



• N 2 c -l 2 ^ 2 9vr 2 1 



(5 j R + J In 2 R + 12 f i? - 1 ) In R 



N 2 s -wp 2 Q 2g 2 

+ 12 (^ + ^) + ( 3i?2 + 2 ( ln (!-^) In 2 - In (i? + 1) In 2 i? 

- 2 Li 2 (-R) InR + 2 Li 2 (i?) In R + 2 Li 3 (-#) - 2 Li 3 (i?))] , (4.10) 

where R = Q1/Q2 ■ When Q 1 = Q 2 , the expression (I4.10p simplifies to 

N 2 -1 9ir 2 
M 00 = is-jjr- a 2 s a em f 2 p ^(14C(3) - 12) . (4.11) 

4.2 Results for differential cross-section 

The formulae for M^; 2 obtained in secH] permit us to evaluate the magnitudes of cross- 
sections (I3.17P of the diffractive double rho production for different helicities of virtual 
photons. In our estimates we use as a strong coupling constant the three-loop running 

MQ1Q2) with A Ks = 305 MeV ( see ' e -g- MS 

In Figj5] we display the t-dependence of the different 7l t 1l,t ~ > Pl Pl differential 
cross-sections for various values of Q = Qi = Q 2 . 

We first note the strong decrease of all the cross-sections when Q\ 2 increase. For LL, 
this follows from an obvious dimensional analysis, since 

(for Qi = Q 2 = Q), in agreement with (14.111) . 

Secondly, all the differential cross-sections which involve at least one transverse pho- 
ton vanish when t = t m i n . It is due to the vanishing of the function Q for r = (see 
(13. 6p ). Physically, this fact is related to the s-channel helicity conservation at t — t min . 
Indeed, since the t-channel gluons carry non-sense polarizations, helicity conservation 
occurs separately in each impact factor. 

In Figl6j we show the shape of the integrands M Al A2 of the various amplitudes -MaiA 2 
as a function of z\ and z 2 , as they appear in formulas (I3.11[ 13.131 and 13.151) : 

Moo = z\ z\ <t>{zi) z 2 z 2 (j)(z 2 ) M 00 (2i, z 2 ) , (4.12) 



Running of as is in principle a subleading effect with respect to our treatment. Nevertheless, nu- 
merically, as we discuss in sec l5.3[ the dependence of our predictions for the rates in e + e~ scattering on 
a choice of a s is negligible at Born order, but is more subtle when LO BFKL corrections are taken into 
account. 
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\t-t min \ {GeV 2 ) 




1 2 3 4 5 

\t-t min \ (GeV 2 ) 

Figure 5: Differential cross-sections for the process 7 £ r Ilt ~ ¥ P^L p\- The s °lid curve 
corresponds to the 7272 mode, the dotted one to the 72 7t mode, the dashed and the dashed- 
dotted ones to the 7^7^/ modes with respectively the same T = T' and different T 7^ T" 
transverse polarizations. The different figures (a), (b), (c) correspond to different values of 

Qi = Qi- 
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Figure 6: Shape of the amplitudes Moo, — M + o, M ++ , — M_| as functions of z\ and z 2 , for 

-t = 0.16 GeV 2 and Q 1 = Q 2 = 1 GeV. 



and for Aj = +, — 




(zi - z x ) (j)(zi) z 2 z 2 4>(z 2 ) M Al0 (zi, z 2 ) , 

(zx - Zi) <f>(zi) (z 2 - z 2 ) <p(z 2 ) M AlA2 (zi, z 2 ) . 



(4.13) 
(4.14) 



Ma 1 a 2 ( 2; 15 #2) is symetric under (zi <-> z~i) for a longitudinal polarization Aj = (cf. 13.41) 
and antisymetric under (z^ <-> Zi) for a transverse polarization Aj = +, — (cf. 13.61) ; thus 
the factors z~i for Aj = and — z~i for Aj = +, — ensure the symmetry of M AlAa under 
(zi <r-> Zi) as we can see on Figj6j Because of the p° L mesons distribution amplitudes 
(f)(zi), Ma 1 a 2 (-2i, z 2 ) vanishes for any polarization in the end-point region. Consequently 
the case of a transverse polarization vanishes in the central region z% = z% — 1/2 and also 
in the end-point region Zi close to or 1, so that it restricts the available Z{ phase-space 
and reduces the resulting differential cross-section, in agreement with the dominance of 
longitudinal photons (helicity conservation) in the process Jlt Ilt ~* Pl Pl- 
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The amplitude involving at least one transverse photon has a maximum at low —t 
value with respect to Q\Q2- The Fig E] corresponds to —t = 0.16 GeV 2 which is a typical 
value for the region where the cross- sections with transverse photons in Fig f5] are maximal. 




Figure 7: Shape of the amplitude -M+_, for -t = 0.01 GeV 2 (left) and -t = 0.8 GeV 2 (right), 
with Q 1= Q 2 = 1 GeV. 



A peculiarly characteristic shape appears in the amplitudes with two transverse pho- 
tons, as shown in bottom panels of FigJBl When the value of t changes towards t m i n 
the peaks become very narrow, as shown in the left panel in FigJ7]for M + _. For t very 
close to t min they are practically concentrated only on the boundary which leads to the 
vanishing of the amplitude. On the other hand, when the value of t increases and leaves 

the maximum of cross-sections the peaks in Fig j6] decrease and spread, as shown for M H 

in the right panel of FigJTl 

In the case of LT polarizations, the shape of the amplitude M + o, which contains only 
one factor [zi — Zi), is shown in the right upper panel of Figj6j Its comparison with the 
upper left panel of Figj6l showing the shape of the M o amplitude, leads to the conclusion 
that M +0 shares some properties with M + _ and M o- In particular, the presence of a 
transverse polarization leads to the vanishing of M +0 at t — t m i n . On the other hand, the 
presence of a longitudinal polarization increases the cross-section at small values of t. As 
a consequence of the competition of these two mechanisms, the maximum of the cross- 
section determined by M +0 is located closer to t min than in the case of the cross- section 
given by M + _. This is illustrated in FigJB] which shows the t-dependence of the various 
differential cross-sections in log-log scale. 

Third, in FigJHJ we display the t-dependence of the ^* LT Ylt ~^ P°l Pl differential 
cross-sections for Q = Q 1 = Q 2 = 1 GeV up to values of —t much larger than photon 
virtualities Qi, where t plays the role of the dominant hard scale in our process. Of course, 
in such a kinematical region the cross-section are strongly suppressed in comparison with 
the small t one. Nevertheless, Figj9] illustrates the expected fact that the hierarchy of 
cross-sections is different in two regions: at large t, the 7^ 7^ — * p° L p° L cross-section 
dominates over the one of 7^ 72 — > p° L p° L , which is the dominant cross-section at small t, 
since the virtual photons are almost on shell with respect to the large scale given by t. 
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0.01 0.1 1 10 

\t - t min \ (GeV 2 ) 

Figure 8: Differential cross-sections for the process T T — > p\ p\ , for small value of t. 
The solid curve corresponds to the 7^72 mode, the dotted one to the 72 7^ mode, the dashed 
and the dashed-dotted ones to the 7^7^; modes with respectively the same T = T' and different 
T ^ T' transverse polarizations, for Qi = Q2 = 1 GeV. 

To conclude this subsection, we note that all the above cross-sections are strongly 
peaked in the forward cone. The phenomenological predictions obtained in the region of 
the forward cone will practically dictate the general trends of the integrated cross-sections. 
This fact is less dangerous than for the real photon case since the virtual photon is not 
in the direction of the beam, and thus the outgoing p mesons can be tagged. The only 
difficulty has to do with the tagging of the outgoing lepton, since the cross-section is 
dominated by small (hard) values of Q\ 2 . In this section we did not modify cross-sections 
by taking into account the virtual photon fluxes, which would amplify both, the dominance 
of small Q 2 region as well as the small yi domain, characteristic for very forward outgoing 
leptons. This is discussed in section [5j In particular, it will be shown that the differential 
cross-sections are experimentaly visible and seems to be sufficient for the t— dependence 
to be measured up to a few GeV 2 . 

Note also that at this level of calculation there is no s-dependence of the cross-section. 
It will appear after taking into account triggering effects and/or BFKL evolution. 

4.3 Quark exchange contribution to the cross-section 

The process (II. 3p described above involves gluon exchanges which dominate at high en- 
ergies. However, at lower energy, the process can be described by double quark ex- 
change. This was investigated in [9], in the case t = t min . Fig{T0] shows in particu- 
lar the diagrams which contributes to the amplitude M# (see FigfS]) for the process 
7l(<7i) 7i(92) —> Piih) Piih) ■ Using the Eqs.(12) and (13) of [H] together with the 
asymptotical p° L distribution amplitude (14.91) one obtains the scattering amplitude for the 
photons longitudinally polarized 
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Figure 9: Differential cross-sections for the process 7£ T ^* LT —> p\ p° L , up to asymptotically 
large t. The solid curve corresponds to the 7£7£ mode, the dotted one to the 7£7^ mode, the 
dashed and the dashed-dotted ones to the 7^7^, modes with respectively the same T = T' and 
different T / T' transverse polarizations, for Q\ = Q2 = 1 GeV. 
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In the large s limit, one respectively gets 
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Figure 10: Feynman diagrams contributing to M# in the case of longitudinally polarized virtual 
photons. 
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Other amplitudes vanish at t = t r 

These expressions should be compared with the corresponding 2 gluons exchange con- 
tributions discussed in the previous sections. The LL amplitude is almost constant around 
t = t min , and given by (14.101). 

The TT amplitude (14.11) behaves as 

M 9 T % ^-ia\s ^1 al a em f p , (4.19) 

where the constant a = 253.5 is extracted from a numerical fit. 

The Eqs. (14.151 - 14. 19j) confirm the well known fact that in the Regge limit the two 
gluon exchange dominates over the double quark exchange. Indeed, the comparison of 
expressions (I4.17P and (I4.18P with formulas (I4.10p and (14.191) shows that gluonic contri- 
butions are proportional to s (in agreement with the usual counting rule < g Scr ^ Af + 1 ) where 
N is the number of t channel exchanged particles of spin <7j). In the case of longitudinally 
polarized photons which does not vanish at t m i n , and for the same photon virtualities 
Q\ = Q\ = Q 2 ■, let us consider the ratio 

M% _ 32 (Ql + gg) Q 2 , . 

LL ~ MH~ 28C(3)-24 sa/ 1 ' 

For a typical value of Q 2 = 1 GeV 2 , as soon as s (~ s 7 * 7 *) is higher than 4 GeV 2 , this 
ratio is bigger than unity, which at first sight seems to be always the case for ILC (14.101) 
would thus completely dominates with respect to ( 14.151) . by several orders of magnitudes. 
In fact, s 7 * 7 * can reach such low value as 4 GeV 2 , because of the outgoing energy carried 
by the outgoing leptons and the strong peak of the Weizsacker- Williams fluxes at small 
7* energies. We discuss this effect in section 15.31 at the level of the e + e~ process, after 
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Figure 11: The LDC project, with the BeamCal forward electromagnetic calorimeter. 



performing the phase-space integration of the differential cross-section at t min . It will be 
shown that nevertheless the quark contribution is really negligible in almost all the ILC 
phase space. 

In the case of the two gluon contribution with transverse virtual photons ( 14.191) which 
vanishes at t = t min , its dominance over the corresponding quark contribution (14.161) 
appears very rapidly when \t — t min \ starts to increase, and persists in the whole essential 
region of the phase space (remember that ( 14.191) is peaked at t — t min = k O.OlGeV 2 where 
k is of order 1-10). This dominance will also be discussed in more detail in section I5TB1 at 
the level of the e + e~ process. 



5 Non-forward Born order cross-section for e + e —> 

e+e ~PL PL 

5.1 Kinematical cuts for the phase-space integration 

Our purpose is now to evaluate the cross-section of the process e + e~ — > e + e~ p\ p\ in the 
planned experimental conditions of the International Linear Collider project [20J. For the 
detector part, we chose to focus on the Large Detector Concept [21], and use the potential 
of the very forward region accessible through the electromagnetic calorimeter BeamCal 
which may be installed around the beampipe at 3.65 m from the interaction point. The 
LDC is illustrated in FigE] 

The cross-section which takes into account all the kinematical constraints, which are 
explained below, is given by 

dQ\ I dQl l dyi dy 2 —— — - — ,(5.1) 



dt J Q 2 . J Q 2 . J e JQ1Q2 dtdy 1 dy 2 dQjdQ% 

^lmzn ^°™"" 
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with Qi m i n = 1 GeV, Qi max = 4 GeV, e = 10 6 and y max = 0.6. 

The cross-section ( 15.11) can be evaluated combining the cross-section formulae (12. 2p . 
(13.171) and the results of section 14.21 for the helicity scattering amplitudes. 

The important feature of the formula (15. ip is that the dominant contribution for the 
e + e~ — > e + e~ p° L p° L process is strongly peaked at low Qi. The integration over Q iy yi is 
peaked in the low yi and Qi phase space region due to the presence in (12.21) of \/{yiQl) 
factors coming from the Weizsacker- Williams fluxes, and thus amplifies this effect. We 
show below that this dominant part of the phase space is accessible experimentaly using 
the BeamCal calorimeter. 

The integration domain in (15. ip is fixed by the following considerations. In the lab- 
oratory frame, which is also the center of mass system (cms) for a linear collider, the 
standard expression for the momentum fractions which respect to the incoming leptons 
and for the virtualities of the bremsstrahlung photons are, respectively, given by 

Vl = % — 1 11 } and Qi = AEE' i sm 2 {9 i /2), (5.2) 

hi 

where E is the energy of the beam, while E[ and 0j are respectively the energy and the 
scattering angle of the out-going leptons. At ILC, the foreseen cm energy is -y/i = 2E = 
500 GeV. The experimental constraint coming from the minimal detection angle 9 m i n 
around the beampipe is given by 9 max = n — 9 min > 9i > 9 m i n and leads to the following 
constraint on y^ 

»>'M = 1 - .JiLm - (5 - 3) 

where the constraint on the upper bound of yi coming from 9 max is completely negligible 
at this cm energy. 

The condition on the energy of the scattered lepton E max > E[ > E min results in 

E m i n E max . » 

Vimax = 1 TT" > Vi > 1 ( 5 - 4 ) 

Moreover we impose that s 7 . 7 * = y^s > cQiQ 2 (where c is an arbitrary constant of the 
order 1) which is required by the Regge kinematics for which the impact representation 
is valid. In subsection 15.31 we show that this constant c can be adjusted to choose bins of 
data for which also in the case of e + e~ scattering the contribution with quark exchanges 
(discussed in sec J4.3p is completely negligible. 

We arbitrarily choose Qi to be bigger than 1 GeV as it provides the hard scale of 
the process which legitimates the use of perturbation theory. Qi max will be fixed to 4 
GeV, since the various amplitudes involved are completely negligible for higher values 
of virtualities Qi values (see section 14.21) . The constraints on yi m in discussed so far are 
summarized by conditions 

yimin = max (f(Qi), 1 - ^p^j and y 2m in = max (f(Q 2 ), 1 - ^jp, ■ ( 5 - 5 ) 



19 







! L 




GeV 


50.S( 
38.14 

25.43 

12.71 

lo 


; 
■ 

L 


5- 

o- 

-5 - 

it) 


— . 














P 






,,,,,,, 
-5 a 

v X, CHL 




S 1 






Figure 12: Energy deposit in GeV in each cell of the BeamCal detector due to beamstrahlung. 

Further simplifications of conditions (15.51) can be done by taking into account that the 
only condition on the maximal value of energy detection of the scattered leptons comes 
from kinematics, i.e. E max = E, and some specific features of the planned detector. 

The BeamCal calorimeter in the very forward region allows in principle to detect parti- 
cles down to 4 mradH More precisely, it measures an energy deposit for an angle between 
4 mrad and 26 mrad. But this detector is also polluted by the photon beamstrahlung, 
specialy for very small angles (see Fig JT2l . We assume a non ambiguous identification for 
particles whose energies are bigger than 100 GeV. More precisely, the efficiency of detec- 
tion of an electron depends on its energy and becomes less ambiguous when the energy 
increases. It is above 70 % in the part of the phase space which dominates the cross-section 
(small yi, corresponding to E\ ~ Ei). A precise evaluation of this efficiency would require 
to set up a Monte Carlo simulation for the beamstrahlung contribution, which is beyond 
the scope of this paper. This sets the maximal value of yi to yimax = 1 — = 0.6 with 
E min = 100 GeV and E = 250 GeV. 

Such a big value of E min can be considered as surprisingly high and could lead to a 
strong reduction of the allowed phase space. In principle one could enlarge the phase 
space by taking into account particles whose energies E[ are between 100 GeV and 20 
GeV with angles 9{ bigger than 10 mrad (see Fig JT2|) . but the contribution of this domain 
is negligible (see FigJT3l) since the lower bound of yi (see Eq. ( 15.51) ) prevents us to reach the 
small values of yi and Qi which give the dominant contribution to the cross-section. We 
safely neglect the contribution of this region of phase space and assume in the following 
Emin = 100 GeV and 9 min = 4 mrad. 

Thus, with 9 min = 4 mrad and y/s = 500 GeV, we have stan 2 (# min /2) = 1 GeV 2 , 
which means that f(Q) < for Q 2 > lGeV 2 . The relations (ET5]) . with E max = E, reduce 
to only one condition y^min — ■ This has to be supplemented numerically with the 

3 Note that in order to get access to any inclusive 4J or diffractive high s processes, one needs very 
small 9 m in- To reach values of 9 m in of a few mrads represents an important technological step which was 
not feasible a few years ago. 
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Figure 13: y\ integration domain for 9 min = 10 mrad, E min = 20 GeV and E max = 100 GeV. 

condition y\ m i n = e, where e is a numerical cut-off: although, because of the Regge limit 
condition, we have y-\ > Sifia > QiminQimm = 6.610~ 6 which thus provides a natural lower 

' t>X S y 2 — Sy2max 1 

cut-off for yi, nevertheless we choose e = 10 -6 so that it is smaller than the smallest 
reachable value of y% but still non zero. This cut-off has no practical effect, except for 
avoiding numerical instabilities in the integration code. 

The above discussion justifies the various cuts in formula (15. ip . 



5.2 Background in the detector 

BeamCal is an electromagnetic calorimeter which cannot distinguish charges of parti- 
cles. Thus, it is important to check that the cross-sections of any other processes which 
could lead to final states which can be misidentified with the final state of the process 
e + e~ — > e + e~p° L p° L are suppressed. Indeed, the final state of the process e + e~ — > 77p° Pl, 
with photons of the same energy deposit in detector as outgoing leptons, cannot be dis- 
tinguished with the final state of e + e~ — > e + e~p° L p° L . 

We shall argue that the process e + e~ — > 77/9° p\ leads to a cross-section which is 
negligible at ILC. Let us first start with the process e + e~ — > p\p\ illustrated in FigfJJJa), 
studied in Ref . [221 [23] . 

Its differential cross-section behaves typically like 



do 
~dt 



oc 



J p 



em J p 

s 2 rrip 



(5.6) 



with the virtualities of the photons propagators equal to m 2 p 



can be found in [23], if one identifies gv-y = fp^p- 

Now, when considering the competitor process e + e~ - 
additional bremsstrahlung photon as in Fig JT^T b). we get 



More accurate expressions 
11P°l Pi; that is adding two 



da' 



+iiplpl 



ll 



do 1 



+e T e p L p L 



dt dyi dy 2 dQ{ dQ\ dt dy\ dy 2 dQ{ dQ\ 



a emQlQt 

aZgS 2 mp 



(5.7) 
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(a) 

Figure 14: Example of Born order diagrams for the process e + e~ 
e + e~ — > 7 7P° Pl process (b). 



2= Pl 



(6) 



Pl Pl ( a ) an d f° r the 



which is suppressed at ILC energies, and would be of comparable order of magnitude only 
for colliders with cm energy of the order of a few GeV. 



5.3 Results for cross-section 

We now display in Fig{15] the cross-sections — -r. as a function of t for the 

different polarizations, which are plotted after integrating the differential cross-section in 
(15. ip over the phase space considered previously. We made the following assumptions: we 
choose the QCD coupling constant to be ^(VQiQz) running at three loops, the parameter 
c = 1 which enters in the Regge limit condition and the cm energy = 500 GeV. Figfl5l 
shows for e + e~ scattering the same differential cross-sections related to different photon 
helicities as Figj5j We see that the shapes of corresponding curves are similar although 
they lead to quite different values of cross-sections. The cross-sections corresponding to 
photons with at least one transverse polarization vanish as in the 7*7* (cf subsection 14.21) 
case at t — t m i n . Similarly, each of them has a maximum in the very small t region. These 
maxima are shown more accurately on the log-log plot in FigJT6l 

At this point one technical remark is in order. By looking into the upper plot in FigfLSl 
related to the .Moo amplitude, one sees that the points corresponding to nonzero \t — t min \ 
approach smoothly the point on the axis \t— t m i n \ = 0. This point \t— t m %n\ = is of special 
interest because it gives the maximum of the total cross-section (since the transverse 
polarization case vanishes at t m i n ) and then practically dictates the trend of the total 
cross-section which is strongly peaked in the forward direction (for the longitudinal case) 
and strongly decreases with t (for all polarizations), as shown already at the level of the 
7*7* cross-sections in subsection 14.21 Due to numerical instabilities, the differential cross- 
section at \t — t m i n \ = must be evaluated in a different way than those for \t — t min \ ^ 0, 
i.e. by the use of expression 1 14. 10p in which the integration over Zi was already done in 
the analytic way. Since Eq. (14.101) involves several polylogarithmic functions its structure 
of cuts is quite inconvenient for further numerical integration over variables yi and Qi. In 
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Figure 15: Cross-sections for e + e~ — > e + e~ p\ p\ process. Starting from above, we display the 
cross-sections corresponding to the 7^7^ mode, to the 727^ modes, to the 7^7^/ modes with 
different T ^T' and finally to the 7^7^/ modes with the same T = T' . 



order to overcome this technical problem it is useful to rewrite (I4.10p by the use of Euler 
identity [24J in the form 



N 2 -l , „ 9 9vr 2 



A^oo = -is c Ar9 a 2 a em f 2 



1 -2 



6 [ R+ —)]n R (5.* 



+ 12 \R--jhiR + 12 [R+rJ + [3R 2 + 2 + — j - Li 2 (1 - R) j \nR 

- In (R + 1) In 2 R - 2 Li 2 (-R) In R + Li 2 (R) In R + 2 Li 3 (-R) - 2 Li 3 (R))] . 

since now the imaginary terms only come from Li 2 (R) and Li 3 (R) along their cuts, which 
cancels among each other analytically. Therefore, one can safely use their real part in a 
numerical fortran code as defined in standard packages. 

The ILC collider is expected to run at a cm nominal energy of 500 GeV, though it 
might be extended in order to cover a range between 200 GeV and 1 TeV. Because of 
this possibility, we below discuss how the change of the energy in cms influences our 
predictions for the cross-sections measured in the same BeamCal detector. Furthermore, 

we discuss the effects of our various assumptions on the cross-section — j e PLPL at 
the point t min , and consequently on the behaviour of the total cross-section. 

FigHT] shows the cross-section at t min as a function of the cm energy ^Js for different 
choices of strong coupling constant a s . To see the sensitivity of our predictions to these 
choices, we plot the cross-section at t min in two cases: the blue curve corresponds to 
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Figure 16: Cross-sections for e + e~ — ► e + e~ p° L p\ process, in log-log scale. Starting from above, 
we display the cross-sections corresponding to the 72 7^ mode, to the 72 7t modes, to the Jt^t" 
modes with different T 7^ T' and finally to the 7^7^; modes with the same T = T'. 

Q-sWQiQi) running at one loop and the red one to asiVQiQz) running at three loops. 
The curves in FigJTT] are very close to each other, which leads to a small uncertainty on 
the total cross-section as we will see in the following. 

The shapes of plots in FigJTT] distinguish two different domains: if the planned cm 
energy range y/s is lower than 500 GeV, the function f(Qi) (cf. equation (15.31) ) appearing 
as a constraint on the minimum value of yi in the phase space integration domain does 
not play any role at Q min = 4 mrad. Thus the cross-section increases with y/s between 200 
and 500 GeV. Because of the condition we assumed on the minimal value of the energies 
of the scattered leptons in the section 15.11 the yi integration domain becomes very narrow 
(cf. equation (15.41) ) when y/s goes to 200 GeV and leads to a strong decreasing of the 
cross-section at this cm energy. Note that if y/s becomes bigger than 500 GeV, f(Qi) 
will cut the small yi region (which contribute mainly because of the Weizsacker- Williams 
photons fluxes) when y/s increases. Thus the cross-section falls down between 500 GeV 
and 1 TeV. This is due to the limitation caused by the minimal detection angle offered by 
the BeamCal calorimeter, which is thus optimal for our process when y/s = 500 GeV. This 
effect on f{Q%) could be compensated if one could increase the value of Qi but this would 
be completely suppressed because of the strong decreasing of the amplitude with Qi. The 
above discussion leads also to the conclusion, that although the Born order cross-sections 
do not depend on s, the triggering effects introduce an s-dependence of the measured 
cross-sections. 

FigiH] shows the cross-section at t m i n for different values of the parameter c which 
enters in the Regge limit condition s 7 . 7 * = y\ yi s > cQ\Q2- The value of the param- 
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Figure 17: Cross-sections for e + e~ — > e + e~ p\ p\ at t = t m in for different a s : the blue and red 
curves for a s running respectively at one and three loops, with c = 1. 

eter c controls the dominance of gluonic contributions to the scattering amplitude: the 
increase of c should lead to suppression of quark exchanges. To see that we display 
the quark contribution in the same bins: we use the usual phase-space for the process 
e + e _ — > e + e~p° L p° L (cf. Eq. fl2.2p ) with the expressions of the amplitudes (14. 15j) and 
(14.161) . and perform their numerical integration on y^, Qi with the same cuts as in the 
two gluon exchange process (cf. section [5TT|) . For each value of c we plot the three curves 
corresponding to the two gluon exchange process and the quark exchange processes with 
longitudinal and transverse virtual photons. 

A technical remark is in order when performing this integration numerically. The 
equation (14.161) is not divergent when Q\ — > s because this limit is only valid if s(l — 
Q\/s){l — Q\j ' s) is finite and positive since this term corresponds in our notation to 
the cm energy of the virtual photons. In order to avoid numerical instabilities we add 
the condition yvg2S > Qf, Q\ to the Regge limit condition. We can check that this 
supplementary constraint does not change our results for the other contributions, namely 
for the two gluon exchange and the quark exchange with longitudinal virtual photons 
processes. 

As expected, the quark contribution is suppressed when increasing c and becomes 
completely negligible as soon as c exceeds 2. All above discussion concerned the case 
t = t m in which determines the general trend of the cross-section in the non forward case. 
Because of that we hope that above conclusions are also valid at the level of the integrated 
over t cross-section. Thus, we omit bellow the quark exchanges. 

We finally obtain the following results for the total cross-section integrated over t. We 
shall show three different predictions which differ by the choice of the definition of the 
coupling constant and by the choice of the value of the parameter c controlling the gluon 
dominance. First we choose ctsiVQiQ?) running at three loops, the constant c = 1, the 
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Figure 18: Cross-sections for e + e~ — ► e + e~ p° L p\ at t = t m i n for different values of the param- 
eter c: the red [black] curves correspond to c = 1, the green [dark grey] curves to c = 2 and and 
the yellow [light grey] curves to c = 3. For each value of c, by decreasing order the curves corre- 
spond to gluon-exchange, quark-exchange with longitudinal virtual photons and quark-exchange 
with transverse virtual photons. 



cm energy y/s = 500 GeV and we obtain (up to numerical uncertainties): 

a LL = 32.4/6 (5.9) 
a LT = 1.5 fb 
a TT = 0.2 fb 

and then 

a™ = 34.1/6. (5.10) 

With a nominal integrated luminosity of 125 fb -1 , this will yield 4.26 10 3 events per year. 

Secondly, with the choice of a s (VQiQ2) running at one loop, the constant c = 1 and 
the cm energy yfs = 500GeV, we obtain: 

a LL = 33.9 fb (5.11) 
a LT = 1.5 fb 
a TT = 0.2 fb 

which leads to 

a Total = 35.6 fb. (5.12) 

As expected, we see that the transition from three to one loop changes very little the total 
cross-section. This result will yield 4.45 10 3 events per year with a nominal integrated 
luminosity of 125 fb -1 . 
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In the third choice, we choose <y s (VQiQ2) running at three loops, the same cm energy 
y/s = 500GeV and the constant c = 2 (for which as previously discussed quark exchanges 
are completely negligible) and we get: 



a LL 


= 28.1 fb 




= 1-3/6 




= 0.2/6 



which imply 

a™ = 29.6/6. (5.14) 

This result will yield 3.7 10 3 events per year with a nominal integrated luminosity of 
1251b" 1 . 

Finally, we also consider the same assumptions as the previous ones except for the 
value of the constant c which is now set to c = 10 in order to consider a more drastic 
Regge limit condition and we obtain: 



which leads to 

Total 



= 19.3/6 (5.15) 
= 0.9/6 
= 0.11/6 



20.3/6. (5.16) 



This result will yield 2.5 10 3 events per year with a nominal integrated luminosity of 
125fb _1 . Thus, this shows that the precise way one implements the restriction of the 
kinematical phase space to the domain of applicability of the impact representation does 
not dramatically change the number of events. 

All the prediction above were obtained using the asymptotical DA. In order to see 
the sensitivity of this assumption on our results, we do also the calculation using the DA 
( 14. 9 p within different models. The choice of the DA of Ref.[25] with a<i = —0.1 and a$ = 
gives 4.2 10 3 events per year, while the choice of the DA of Ref.[26] with 02 = 0.05 and 
a 4 = gives 4.3 10 3 events per year. 

In summary of this part we see that our predictions are quite stable when changing 
the main parameters characterizing the theoretical uncertainties of our approach. 

The obvious question which appears now is how our predictions summarized by Eqs. 
(15. 10115. 141) will change by the inclusion of the BFKL resummation effects. Generaly BFKL 
evolution increases strongly values of cross-sections, which means that usually the results 
obtained at Born approximation can be considered as a lower limit of cross-sections for 
p-mesons pairs production with complete BFKL evolution taken into account. Although 
the complete analysis of BFKL evolution for our process is beyond the scope of the present 
paper, we would like to finish this section with a few remarks on possible effects caused 
by the BFKL evolution. 
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We consider below only the point t = t min and we restrict ourselves to the leading 
order BFKL evolution^. Of course such an estimate should be taken with great caution 
since it is well known that LO BFKL overestimates the magnitude of corrections. 
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Figure 19: Cross-sections for e + e~ — > e + e~p° L p° L with LO BFKL evolution at t = t m i n for 
different a s : the upper and lower red curves for a s running respectively at one and three loops 
and the green one (the middle curve) for a s = 0.46. 



In the Fig{19] we show the corresponding cross-section at t min as a function of y/s, for 
different choices of a s \ we considered a s running at one and three loops (red curves) as 
in the previous discussion for the two gluon exchange and we also used a fixed value of 
a s (green curve) corresponding to the three loops running coupling constant at a typical 
virtuality Q — 1.1 GeV. We have used the expression of the BFKL amplitude [8] for the 
forward case in the saddle point approximation, namely 

Af + + n ■ 5 ^ 9(N?-l)a 2 s aemf 2 p e 41n2 ^ y / In 2 R \ 

A(s, t = t mm , Ql , Q 2 ) ~ ls « ^ — m ^^y^ {-Y^CW) ' 

(5.17) 

with the rapidity Y = ln(^f^), a s = ^aJv^C}^) and R=&. The plots in FigHSI 
are obtained by assuming that the constant d in Eq. fl5.17l) . which at LO is arbitrary and 
of order 1, is chosen to be 1. The factor exp(41n2 a s Y) explains the enhancement of the 
sensitivity to the choice of a s compared to the one in the Born two gluon exchange case, 
since 4 In 2 Y takes big values for ILC rapidities Y. For the same reasons as discussed 
earlier in this section, the function f(Qi) does not appear for yfs lower than 500 GeV; 
the LO BFKL cross-section then grows exponentially with s in this domain. The effect 
of f(Q) starting from 500 GeV gives an inflexion point of the curves and a maximum 
beyond 500 GeV; then the curves decrease until ITeV. 



4 Note that related studies with a hard scale provided by t and not by Q 2 were performed in [27] 
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200 400 600 800 1000 

(GeV) 

Figure 20: LO BFKL cross-section for e + e~ — ► e + e~p° L p° L at i = i m j„ for different values of 
the parameter c': by decreasing order, the curves correspond to c' = 2, c' = 1 and d = 0.5 . c is 
fixed to be equal to 1. 



The effect of varying the parameter d in the BFKL prediction is illustrated in Fig]20l 
As expected, it has a strong effect in the order of magnitude of the differential cross- 
section, since the rapidity is very high and thus leads to a large value of the factor 
exp(41n2 a s Y), which is highly sensitive to the precise definition of the rapidity. 

Comparing the order of magnitude of Born cross-section (Fig{T7] and [T8|) with cross- 
sections provided by the LO BFKL evolution (Fig {19] and Fig j20|) . one could be astonished 
by the fact that they differ by several orders of magnitudes. From previous studies at the 
level of 7*7* [5]. [15]. the NLO contribution is known to be between LO and Born order 
cross-section. Thus, at the level of the e + e~ process, such a large magnitude for the LO 
BFKL cross-section will be suppressed at NLO, leading to a more realistic estimate. 

The above discussion about BFKL enhancement was restricted to the forward case 
t = tmin- in the non-forward case, the phase space region with small t values dominates 
the cross-sections. The obtained hierarchy between cross sections in Born approximation 
for different photon polarizations will presumably still be valid when including BFKL 
evolution at any order of resummation (LO, NLO, etc.). Indeed the argument given in 
section H~2l for Born order and on which this hierarchy is based, only relies on the s-channel 
helicity conservation. Technically, it is based on the impact representation which is valid 
beyond Born and/or LO approximation. 

The comparison of Figs nTlfTSI with Figs fTTflliZUl leads to the conclusions that the BFKL 
evolution changes the shape of the cross-section: when increasing y/s from 500 GeV 
to 1 TeV, the two gluon exchange cross-section will fall down, while with the BFKL 
resummation effects, the cross-section should more or less stay stable, with a high number 
of events to be still observed for these cm energies. 
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6 Conclusion 



The present study should be considered as a continuation of our previous investigations 
[6],[Zj,[8] for the production of two p° -mesons in the scattering of two longitudinally 
polarized virtual photons. The diffract ive production of a meson pair is one of the gold 
plated processes which permit clean studies of the BFKL dynamics at ILC Our main 
motivation in the present work was to estimate, in the Born approximation, the cross- 
section for production of p^-meson pairs in the e + e~ collisions occurring in the kinematical 
conditions of future ILC. For this aim, we first calculated contributions, missing up to 
now, which involve the helicity amplitudes with transversally polarized virtual photons. 
This was done in a mostly analytic way, by the use of techniques developped in Ref. [Tj. 
Having done, we calculated the cross-section for the electroproduction of p°-meson pairs 
which takes into account kinematical cuts imposed by the LDC design project for the 
BeamCal detector. By assuming a nominal value of the integrated luminocity, we predict 
(in the numerical analysis of cross-sections) a production of at least 4 10 3 meson pairs per 
year, a value which is sufficiently large to ensure a reliable data analysis. 

We discussed a possible background process in the BeamCal detector which can iden- 
tify in a misleading way an outgoing lepton with a photon. We predict that the cross- 
section for such a background process is negligibly small at ILC energies. 

Finally we discussed theoretical uncertainties of obtained estimates. There are two 
main sources of them. The first one is related to the assumptions we have made to 
characterize the Regge limit and the particular role played by the parameter c; we also 
observe a sensitivity of our results on the choice of the running coupling constant. 

The second source of theoretical uncertainties of our estimates is related to taking into 
account effects of the BFKL evolution. Generally, an inclusion of the BFKL evolution 
increases significantly the cross-section, as one sees from comparison of Fig {T7| with Fig {19] 
obtained within the LO BFKL approach. On the other hand it is known that this increase 
of predictions for cross-sections is smaller if the BFKL evolution is considered at the 
next-to-leading order. Because of that we can safely say that our predictions should be 
considered as a lower limit of predictions which are obtained by taking into account BFKL 
effects at NLO. We hope to consider this issue in our future publications. 

In principle, the same techniques can be applied for the description of processes in- 
volving other final states, both with positive charge parity exchange in t— channel, e.g. 
J/ty pairs, as well as negative charge parity, e.g. 7*7* — > rj c r} c [2"g] . 
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7 Appendix 



In this appendix we collect the analytical expressions of the coefficients a(r; Q±, Q 2 ; Zi, z 2 ), 
b(r; Qi, Q 2 ; z±, z 2 ) and /(r; Qi, Q2', z\, Z2), as well as all the generic integrals which appear 
in the computation of the Born amplitude in section HI The coefficients a, b and / can 
be expressed as combinations of several generic integrals 

a(r; Qi, Q 2 ; z%, z 2 ) = - [hfn^izi, z 2 ) + l3 MlAta (zi, z 2 ) — h[i, ltH (zi, z 2 ) — l^^^i, z 2 )\ 
^2 

~~2 IAmi/mC 2 !' Z V ~ ^4 WA t 2 (^l) z 2)) 



b{r;Qi,Q2;z 1 ,z 2 ) 



ziz 2 



Z\Z 2 



[z 2 r 2 + n\)(zlr 2 + (2; 2 r 2 + /xf)(^r 2 + /i 2 ) 



2l2 2 



+ 



212:2 



+ 

+ 
+ 



zfr 2 + (z 2 r 2 + /x 2 ,) (zfr 2 + fx 2 ) (z 2 r 2 + /i 2 ) 
z\ Z\ 



r 2 h 



zfr 2 + n\ zfr 2 + \x\ 
: + 2z x z 2 - (zi + z 2 ) 



[h^{z 2 ) - I 2 ^(z 2 ) + r 2 (z 2 - z 2 )h^{z 2 )] + (1 <-> 2)j 



-- + (21 + 22) 



21 - 2 2 



(2l, Z 2 ) + 

(2l, Z 2 ) + 



[-^1^2(21, 2 2 ) + (2l, 2 2 )] 

3 



(21 + 2 2 ) 



(2l, 2 2 ) 



2 2 - 2l 



(21, 2 2 )] 



"^[^MlA^ 1 ' + Ju lf i 2 (zi, 2 2 ) - J"u 1 „„(2l, Z 2 ) - ^„ 1 u 2 (2l, 2 2 )] , 



r 



(7.2) 



f(r;Q 1 ,Q2;z 1 ,z 2 ) 
+ 2 * 2 



+ 



z 2 r 2 + /x 2 z 2 r 2 + /i 2 



2 2 



22 



2|r 2 + /i 2 z 2 r 2 + \i\ 



Z-2 



4 m (2l) + {Z 2 ~ 2 2 ) 



+ 



z 2 r 2 + nl zfr 2 + /i 2 



^2(22) 



+ (2 2 — 2 2 ) [/4 MlM2 (2i, 2 2 )+/ 4MlAt2 (Zi, 2 2 )] + - 



+ 



2 2 z: 2 + n\ zfr 2 + /if 



[^ 2 (2 2 )-/ 2At2 (2 2 )] 



H — n [^1/^2(21, 2 2 ) — h^ lfJl2 (zi, z 2 ) + I^^zi, z 2 ) — I^^zi, z 2 )\ 



(7.3) 



part of which are divergent integrals 



d d k 



(7.4) 
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and some of them are finite 
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J^n^yZii z 2 ) 



d d k 
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(7.15) 



,2 <-> Z) 



,((k - r_ Zl ) 2 + //?)((& - rz 2 ) 2 + « 2 ) {z\r 2 + /if)(^ 2 2 r 2 + /i 2 ) 

Although divergent integrals (17.4117. 8p enter formulas (17. If7.2l and 17.31) . the coefficients 
a fc Qi, Q2; zi, z 2 ), b(r; Q 1 , Q 2 ; Zi, z 2 ) and /(r; Qi, Q 2 ; Zi, z 2 ) are finite expressions. This 
requires a tracing of the cancellation of divergencies within each of these coefficients which 
in principle can be done within a dimensional regularization scheme. Unfortunately this 
procedure is not always practically feasible, especially for the integrals I^m, h mi m 2 or 
/4 mim2 , which apart of massless propagators contain also massive propagators. These 
integrals can be rewritten in terms of simpler, although divergent integrals I 2 and I 2m as 



Ol, Z 2 ) 



1 



2 \ z(r 2 + m\ z 2 r 2 + mi 



^2mim 2 ( z li Z 2 ) 



(7.16) 
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2l 



2 2 



v : 2 r 2 • ///- z 2 r 2 + m 2 1 



(zi, z 2 ) + 



2(z 2 r 2 + m 2 . 



2\ J2m i 



(*0 + 



2(zfr 2 + m\ 



2 \ J 2m 2 



(*0 



'4mim2 



(^1, 2 2 ) 



1 



+ 



2 \(zfr 2 + m\){z 2 z 2 + m|) (zfr 2 + ml)(z 2 r 2 + m 2 ." 
1 



(Zl, Z 2 ) 



2 L 2 r 2 + m 2 + + m 2 ) h + 2 J3m(z) ' 



(7.17) 
(7.18) 



Divergent integrals like I2 and l^m or finite ones like J™ im2 and J2 mi m 2 are directly com- 
puted by standard Feynman techniques. 



The finite integrals J^ mi m 2 an d Jz m are calculated by the use of the change of vari- 
ables defined by the inverse conformal transformation applied to the momenta and other 
dimensional parameters of integrals, see [5] for details. Such a change of variables applied 
to the two-dimensional, UV and IR finite integrals results in a reduction of number of 
massless propagators, a well known property from the studies of conformal field theories 
in coordinate space [29J. 

As a result, we get the following expressions for the coeffcients a(r; Qi, Q 2 ; zi, Z2), 
b(n Qi, Q2] Zi, z 2 ) and f(r; Qi, Q2; z%, z%) which involve only explicitely convergent inte- 
grals 
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The separation of divergent and finite parts in all integrals involved was done within 
the dimensional regularization with d = 2 + 2e. Thus for the simplest integral I 2m the 
result reads 



i2m(o) 



7T 



e(a 2 r 2 + m 2 ) 
or for its massless case 



(1 + e(ln7r - (1) - lnm 2 + 2 ln(aV + m 2 ))) 



2ir 
r 2 e 



(l + e(ln(vrr 2 )-^(l))) 



(7.22) 



(7.23) 



We will use the finite part of i 2 and I 2m to compute the integrals in which they are 
involved. For example we can write the finite integrals K marrih , L maTrib and M maTnb in 
term of the finite part of these integrals since the divergent parts ultimately cancel when 
we express them like 



K, 
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(7.25) 
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The integrals which involve only two massive propagators are finite and are calculated in 
a standard way using Feynman parameters. We obtain 



2m a m b 



d d k 



((k ~ L a ) 2 + m a)((k — zt>y 



mi 



it ^ r 2 (a — b) 2 + m 2 a + ml + \/~X 



(7.27) 



y/X r 2 (a — b) 2 + ml + m 2 — \/X 
where we introduce the notation 

X(x, y, z)=x 2 + y 2 + z 2 - 2xy - 2xz - 2yz , (7.28) 
which enables us to define, for the purpose of our computation, 

A = \(-r 2 (a - 6) 2 , a 2 , (3 2 ) = (a 2 - (3 2 ) 2 + 2(a 2 + (3 2 )r 2 (a - b) 2 + r 4 (a - bf . (7.29) 
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(7.30) 



The In | terms which apparently give divergencies when a or b goes to zero are spuri- 
ous: indeed they cancel because of the symmetrical way the function J ma m b appears in 
afc Qi, Q2; zi, z 2 ), b(r; Qx, Q 2 ; Zi, z 2 ) or /(r; Qx, Q 2 ; z x , z 2 ). 

Let us now consider the J^ m integral. After performing a special conformal transfor- 
mation, namely an inversion on momentum integration variables and other dimensional 
vectors and parameters, a translation and again an inversion, we arrive to integrals with 
smaller number of propagators which are calculated in the standard way. The final result 
reads 
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. (7.31) 



With the same but more tricky approach we can compute the J^ mm integral after per- 
forming a special conformal transformation; we refer to [7] for the complete calculation 
and final expression of this integral (see (A. 66) and (A. 67) of [7]). Starting from this 
result, we now derive another expression in such way that some spurious divergent terms 
(appearing when r 2 = Q 2 , which corresponds to a = or b = 0) explicitely cancel. Thus 
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it allows us to use it in our numerical integration code (cf. sections |4] and [5]). This makes 
the explicit expression (17.331) different from the one given in (A. 67) of [7j. The resulting 
expression for is 
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where appropriate additional In r 2 terms have been introduced in order to write the final 
result as made of logarithms of dimensionless quantities. The difference between these 
two expressions only involves terms proportional to lnr 2 , in accordance to the dimensional 
regularization, since J^p is a divergent integral. However, at the level of the final result 
for J^ x n 2 i which is both UV and IR finite, these additional terms of course cancel each 
other. 
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